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ABSTRACT 



> ; 

V) , Aims. In order to fully understand the physical processes in the magnetospheres of the Magnetic Chemically Peculiar stars, we performed 
OS multi-frequency radio observations of CU Virginis. The radio emission of this kind of stars arises from the interaction between energetic 
electrons and magnetic field. Our analysis is used to test the physical scenario proposed for the radio emission from the MCP stars and to 
derive quantitative information about physical parameters not directly observable. 

Methods. The radio data were acquired with the VLA and cover the whole rotational period of CU Virginis. For each observed frequency 
\ the radio light curves of the total flux density and fraction of circular polarization were fitted using a three-dimensional MCP magnetospheric 
"-■■«». model simulating the stellar radio emission as a function of the magnetospheric physical parameters. 

' Results. The observations show a clear correlation between the radio emission and the orientation of the magnetosphere of this oblique 
rotator. Radio emission is explained as the result of the acceleration of the wind particles in the current sheets just beyond the Alfven 
Q ■ radius, that eventually return toward the star following the magnetic field and emitting radiation by gyrosyncrotron mechanism. The 
' accelerated electrons we probed with our simulations have a hard energetic spectrum (N(E) cc E~ 2 ) and the acceleration process has an 
C/3 , efficiency of about 1CT 3 . The Alfven radius we determined is in the range of 12 - 17 R, and, for a dipolar field of 3000 Gauss at the 
magnetic pole of the star, we determine a mass loss from the star of about 10 12 Mq yr '.In the inner magnetosphere, inside the Alfven 
\ radius, the confined stellar wind accumulates and reaches temperatures in the range of 10 5 - 10 6 K, and a detectable X-ray emission is expected. 
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1. Introduction 

Magnetic Chemically Peculiar (MCP) stars are characterized 
by periodic variability of the effective magnetic field through 
the stellar rotational period. The magnetic field topology of this 
kind of star has been commonly regarded as a magnetic dipole 
tilted with respect to the rotation axis (Babcock 1949, Stibbs 
rT950l . even if MCP stars with a multipolar magnetic field are 
also known (Landstreet l 19901 Mathvs fTWTT l. 

In some cases, MCP stars are characterized by an 
anisotropic stellar wind (Shore et al. 119871 Shore & Brown 
I1990> . as a consequence of the wind interaction with the dipo- 
lar magnetic field (Shore fT987> . The resulting magnetosphere 
is structured in a "wind zone" at high magnetic latitudes, where 
the gas flows draw the magnetic field out to open structures, and 
a "dead zone", where the magnetic field topology is closed and 
the gas trapped. 



Currently, non-thermal radio emission is observed from 
about 25% of MCP stars and the detection rate seems to be 
directly correlated to the effective stellar temperature (Drake 
et al. [T9871 Linsky et al. [T992I Leone et al. [T994T l. Since 
only MCP stars with high photospheric temperature develop 
a radiatively-driven stellar wind, this result reveals the impor- 
tance of a stellar wind in driving the radio emission. To explain 
the origin of the radio emission from pre-main sequence mag- 
netic B stars and from MCP stars, Andre et al. ( 1988 ) proposed 
a model characterized by the interaction between the wind and 
the magnetic field. The radio emission is understood in terms 
of gyrosynchrotron emission from mildly relativistic electrons, 
which are accelerated in "current sheets" formed where the gas 
flow brakes the magnetic field lines, close to the magnetic equa- 
tor plane, and propagate along the magnetic field lines back 
to the inner magnetospheric regions (Havnes & Goertz 1984; 
Usov & Melrose fT992l . 
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According to the oblique dipole model, Leone ([1991 ) sug- 
gested that the radio emission of MCP stars should be also pe- 
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Table 1. Results of the observations. For each observing scan, the mean UT time, the total flux density (Stokes/) at the three 
observed frequency (7 5 , 7g 4 and 7 15 ), the corresponding errors (<x) and the circular polarized flux density, Stokes V (Vg, V% .4 and 
V15), are reported: 



<UT> 


/ 5 (cr) 


V 5 (cr) 


<UT> 


44 CO") 


V SA (cr) 


<UT> 


/is (cr) 


V l5 (cr) 




[mJy] 


[mJy] 






[mJy] 


[mJy] 




[mJy] 


[mJy] 


Date: 1998 June 02 


00:30:45 


2.78 (0.05) 


< 3a- (0.05) 


01:37:00 


3.18 (0.04) 


< 3a- (0.04) 


02:51:50 


4.1 (0.1) 


< 3a- (0.1) 


01:23:05 


3.27 (0.05) 


0.32 (0.05) 


02:36:55 


3.87 (0.04) 


< 3a- (0.04) 


03:51:35 


4.3 (0.1) 


< 3a- (0.1) 


02:22:55 


3.64 (0.04) 


< 3a- (0.04) 


03:38:50 


4.40 (0.04) 


< 3a- (0.04) 


04:51:25 


3.4 (0.1) 


< 3a- (0.1) 


03:22:45 


4.13 (0.04) 


< 3a- (0.04) 


04:36:35 


4.42 (0.03) 


< 3a- (0.03) 


05:51:20 


3.3 (0.1) 


< 3a- (0.1) 


04:22:35 


4.44 (0.04) 


< 3a- (0.04) 


05:36:30 


3.69 (0.03) 


< 3a- (0.03) 


06:51:15 


2.1 (0.1) 


< 3a- (0.1) 


05:22:25 


3.48 (0.04) 


< 3a- (0.04) 


06:36:20 


2.85 (0.03) 


< 3a- (0.03) 


07:44:20 


3.0 (0.1) 


< 3a- (0.1) 


06:22:15 


2.78 (0.04) 


< 3a- (0.04) 


07:30:25 


2.85 (0.04) 


0.22 (0.04) 








07:18:55 


3.96 (0.05) 


0.58 (0.05) 
















08:13:45 


3.57 (0.05) 


0.27 (0.05) 
















Date: 1998 June 07 


02:10:15 


4.87 (0.04) 


0.43 (0.04) 


02:23 


40 


4.85 (0.03) 


0.77 (0.03) 


02:38:10 


3.3 (0.1) 


0.4 (0.1) 


03:08:35 


4.87 (0.04) 


0.37 (0.04) 


03:21 


35 


4.97 (0.03) 


0.56 (0.03) 


03:35:55 


3.0 (0.1) 


< 3a- (0.1) 


04:06:25 


4.68 (0.04) 


0.34 (0.04) 


04:19 


25 


4.52 (0.03) 


0.57 (0.03) 


04:33:45 


3.6 (0.15) 


< 3a- (0.15) 


05:04:50 


3.73 (0.04) 


< 3a- (0.04) 


05:18 


20 


3.50 (0.04) 


0.28 (0.04) 


05:32:45 


3.3 (0.15) 


< 3a- (0.15) 


06:03:35 


3.79 (0.04) 


< 3a- (0.04) 


06:17 


05 


3.88 (0.04) 


< 3a- (0.04) 


06:31:35 


3.4 (0.15) 


< 3a- (0.15) 


07:02:30 


3.98 (0.04) 


< 3a- (0.04) 


07:15 


15 


4.01 (0.04) 


< 3a- (0.04) 








07:54:30 


4.38 (0.05) 


< 3a- (0.05) 
















Date: 1998 June 12 


01:21:35 


3.95 (0.05) 


< 3a- (0.05) 


01:35:05 


4.31 (0.04) 


-0.14(0.04) 


02:47:25 


3.3 (0.1) 


< 3a- (0.1) 


02:19:55 


4.10 (0.04) 


< 3a- (0.04) 


02:32:55 


4.18 (0.04) 


-0.21 (0.04) 


03:45:15 


2.5 (0.1) 


< 3a- (0.1) 


03:17:45 


3.47 (0.04) 


< 3a- (0.04) 


03:30:45 


3.54 (0.03) 


< 3a- (0.03) 


04:44:10 


2.7 (0.1) 


< 3a- (0.1) 


04:16:05 


3.12 (0.04) 


< 3a- (0.04) 


04:29:35 


3.05 (0.03) 


0.16(0.03) 


05:43:00 


3.3 (0.1) 


< 3a- (0.1) 


05:15:05 


3.28 (0.04) 


0.50 (0.04) 


05:28:25 


3.58 (0.03) 


0.39 (0.03) 


06:41:50 


3.7 (0.1) 


< 3a- (0.1) 


06:13:55 


4.38 (0.04) 


0.44 (0.04) 


06:27:20 


4.33 (0.03) 


0.53 (0.03) 








07:09:30 


4.69 (0.04) 


0.43 (0.04) 


07:20:55 


4.67 (0.04) 


0.93 (0.04) 









riodically variable. The confirmation of this hypothesis comes 
from the observations of HD37017 and HD374 79 (Leone & 
Umana fT993l and of HD133880 (Lim et al. 1996 1, performed at 
5 GHz, that revealed the presence of rotational flux modulation. 
In particular, the observed flux light curves of HD37017 and 
HD37479, which are characterized by a simple dipolar mag- 
netic field, show a radio emission minimum coincident with the 
zero of the longitudinal magnetic field, whereas the emission is 
maximum at the magnetic field extrema. 

Among radio MCP stars, only CU Virginis (HD 124224) is 
characterized by a strong 1.4 GHz flux enhancement around 
phases 0.4 and 0.8, with a right hand circular polarization of 
almost 100% (Trigilio et al. 2000] TLL00). This discovery has 
been explained in terms of Electron Cyclotron Maser Emission 
(ECME) from energetic electrons accelerated in the current 
sheets near the Alfven surface (where the thermal plus kinetic 
plasma pressure equals the magnetic one) and reflected out- 
ward by magnetic mirroring in the inner magnetospheric lay- 
ers. 

In a previous paper (Trigilio et al. 2004 TLU04), we devel- 
oped a three-dimensional numerical model of the radio emis- 
sion from MCP stars with a dipolar field. On the basis of suc- 
cessfull simulated 5 GHz flux light curves of HD37017 and 
HD37479, whose study allowed us to confirm the qualitative 



scenario proposed to explain the origin of non-thermal radio 
emitting electrons. We showed that the combined effect of a 
tilted dipolar magnetic topology and the presence of a thermal 
plasma trapped in the "dead-zone" can produce the observed 
phase modulation of the stellar radio emission. 

In this paper we present the multifrequency microwave ob- 
servations of CU Virginis, a well studied Si-group MCP star. 
To investigate the physical conditions of the stellar magneto- 
sphere, the radio light curves of CU Vir are compared with the 
results of the numerical simulations performed using our 3D 
model. We present here an implementation of our code to also 
treat circular polarization. 

2. Observations and data reduction 

The observations were carried out using the VLA 1 in three dif- 
ferent runs, on June 02, 07 and 12, 1998. Each run was ap- 
proximately eight hours long. We observed at four frequencies 
(1.4, 5, 8.4 and 15 GHz) using two independent and contigu- 
ous bands of 50 MHz in Right and Left Circular Polarizations 

1 The Very Large Array is a facility of the National Radio 
Astronomy Observatory which is operated by Associated Universities, 
Inc. under cooperative agreement with the National Science 
Foundation. 
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Fig. 1. Top panel: effective magnetic field of CU Vir versus 
rotational phase (triangles: data by Borra & Landstreet 1980 
squares: data by Pyper et al. 1998). Other panels: phase mod- 
ulation of Stokes / and Stokes V flux density at 15, 8.4 and 5 
GHz. For comparison, the 1 .4 GHz Stokes V (lowest panel) by 
TLL00 is reported. The three days of observations are distin- 
guished with different symbols to point out the repeatability of 
the radio emission in time. 



(RCP and LCP). The observations were performed using the 
entire array (configuration array moving A — > B) alternating all 
the frequencies. 

Because of the w Q.52 d rotational period and with the aim 
to obtain a complete rotational phase coverage, we assumed at 
any frequency an observing cycle consisting of a: 10 minutes 
on source, preceded and followed by «2 minutes on the phase 
calibrator (1354 - 021). To avoid redundancy in the phase cov- 



erage for each frequency, the observing sequence was changed 
in each run. 

The flux density scale was determined by daily observa- 
tions of 3C286. The data were calibrated and mapped using 
the standard procedures of the Astronomical Image Processing 
System (AIPS). The flux density and fraction of circular polar- 
ization were obtained by fitting a two-dimensional Gaussian, 
task JMFIT, at the source position in the cleaned maps (sepa- 
rately for the Stokes / and Stokes V 2 ) integrated over the du- 
ration of each individual scan. As the uncertainty in the flux 
density measurements we assume the r.m.s. of the map. Results 
are summarized in Tabled The measurements performed at 1 .4 
GHz have been presented elsewhere (TLL00). 



3. The radio behavior of CU Vir 

In Fig.[0 for each observed frequency, we report the radio flux 
density measured in both Stokes / and V versus the rotational 
phase. This was computed adopting the ephemeris referred to 
light minimum reported by Pyper et al. J1998l >: 



HJD = 2435 178.6417 + 



d .52067780£ JD < 2446 000 
d .52070308£ JD > 2446 000 



In the lowest panel of Fig. [2 the measures performed in 
Stokes V at 1.4 GHz from TLL00 are also displayed. 

The top panel of Fig.^shows the average longitudinal com- 
ponent of the magnetic field as a function of the rotational 
phase (data are from Borra & Landstreet 1980 and Pyper et 
al. ri99k 

The flux variability with the stellar rotational period sug- 
gests that the radio emission at 15, 8.4 and 5 GHz arises from 
a stable co-rotating magnetosphere and that the radio source is 
optically thick. For Stokes / the radio emission minimum ap- 
pears to be coincident with the null magnetic field (<J> =s 0.4 and 
O « 0.8), whereas the two maxima, clearly visible at 5 and 8.4 
GHz, are related to stellar orientations coinciding with the two 
extrema of the stellar magnetic field (<£ m 0.1 and <t> » 0.6). 
Stokes V is detected at all the observed frequencies close to 
<t> » 0.1, whereas at <t> * 0.6 the circular polarization is de- 
tected only at 8.4 GHz. The circularly polarized flux density 
reaches its positive (RCP) maximum value when the effective 
magnetic field strength is also maximum. This means that ra- 
dio flux emitted in a direction close to the magnetic north pole, 
where magnetic field lines are mainly radially oriented (relative 
to the outside stellar surface), is partially right hand polarized. 
Negative values of Stokes V (LCP) have been detected only at 
8.4 GHz around O * 0.6, when the magnetic south pole is clos- 
est to the line of sight. 

The phase modulation of the Stokes / flux density has been 
observed in other MCP stars at 5 GHz (Leone [T99U Leone 
& Umana fT993l Lim et al.[T996 1. Our dataset represents an ex- 
tension of previous studies for the (almost simultaneous) multi- 
frequency observation of both total and circular polarized flux 



2 VLA measurements of wave polarization state are in agreement 
with the IAU and IEEE orientation/sign convention: positive V corre- 
sponds to right circular polarization. 
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/ 

Fig. 2. Schematic view of the magnetosphere of an MCP star 
assuming a central dipolar magnetic field. The dashed line 
shows the Alfven surface: the largest magnetic field line de- 
limiting the middle magnetosphere is drawn with a dotted line 
outside the Alfven surface. The length / represents the current 
sheets region where the non-thermal electrons are efficiently 
accelerated. The inner magnetosphere is filled with hot thermal 
plasma and cold material accumulates in the form of a torus 
near the stellar surface. 

with a detailed phase coverage of the rotational period. This al- 
lows us to monitor the variations of an MCP star radio spectrum 
through its entire stellar rotation. 

On average the observed radio spectra are flat with a 
slightly negative power law index a « -0.1 (I v oc v"). This, 
together with the low fraction of circular polarization (about 
10%), supports the idea that non-thermal gyrosynchrotron from 
a non-homogeneous source is responsible for the observed ra- 
dio emission. 

4. The model 

Our 3-D model samples the magnetosphere of the star in a three 
dimensional grid, assumes a set of physical parameters, com- 
putes emission and absorption coefficients for gyrosynchrotron 
emission, and solves the equation of transfer for polarized radi- 
ation in terms of Stokes / and V (appendix|A}. From the best fit 
of the observed fluxes, we determine the most probable physics 
of the magnetosphere. 

4.1. Basic picture 

The model was developed (TLU04) on the basis of the physi- 
cal scenario proposed by Andre et al. ( 1988) and Linsky et al. 
( 1992 1. The model involves the interaction between the dipolar 
magnetic field and the stellar wind resulting in the formation 
of "current sheets" outside the equatorial region of the Alfven 



Table 2. Adopted stellar parameters for CU Virginis. 



Distance [pc] 


D 


80 T 


Stellar Radius [R©] 


R 


2.2 a 


Polar Magnetic Field [Gauss] 


Bp 


3000 d 


Rotation Axis Inclination [degree] 


i 


43 d 


Magnetic Axis Obliquity [degree] 


P 


74 d 



ESA fT997l 
a North fT998i 
b Kuschnig et al. fP%9l 
c Pyper et al. (19981 
d Trigilio et al. i2000l 



surface (Havnes & Goertz 1984 >, where the wind particles can 
be accelerated up to relativistic energy (Usov & Melrose 1992 1, 
radiating at radio wavelengths by the gyrosynchrotron mecha- 
nism. 

The adopted geometry is shown in Fig.|2] The layer where 
relativistic electrons propagate is delimited by two magnetic 
lines: the first is the largest closed line at a tangent to the Alfven 
surface in the plane of the magnetic equator and the latter is 
defined by the length of the current sheets where electrons are 
efficiently accelerated. In each point of this shell, the "middle" 
magnetosphere in analogy with Andre et al. 0988), it is as- 
sumed that emitting electrons have a power law energy spec- 
trum, isotropic pitch angle distribution and homogeneous spa- 
tial distribution due to magnetic mirroring. 

Following the magnetically confined wind shock (MCWS) 
model by Babel & Montmerle (1997), the "inner" magneto- 
sphere ("dead zone") is filled by a trapped thermal plasma 
heated to a temperature higher then the photospheric value 
(r phot = 12500 K, Kuschnig et al. [19991 This hot thermal 
plasma progressively pushes the shock front to high magnetic 
latitudes, until its thermal pressure equals the wind ram pres- 
sure. 

The net effect of this mechanism is the filling of the inner 
magnetosphere with hot thermal plasma. If the magnetosphere 
rotates, the density decreases linearly outward, while the tem- 
perature increases linearly (Babel & Montmerle. 11997b . 

To balance the continuous supply of hot matter from the 
polar regions, a cool plasma down-flow is expected from the 
inner magnetosphere toward the stellar surface. This process 
causes the accumulation of cold dense material, which has been 
assumed for simplicity as a homogeneous torus with a circular 
cross-section surrounding the star (Babel & Montmerle 1 1 9971 
Smith & Groote l200Tl . 

4.2. The parameters 

In Table|5]we list the assumed stellar parameters. The depen- 
dence of our simulations on the polar magnetic field B p , incli- 
nation ; and obliquity /3 that are not known with great accuracy 
will be discussed later in Sect. 15.21 

The free parameters of the model are: 

r\ average radius of the Alfven surface in R* units; 
/ thickness of the middle-magnetosphere at r& in R* units; 
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Symbol 


Range 


Simulation step 


Free parameters - step 1 


Size of Alfven radius [R„] 


r A 


6- 


-30 


Ar A 


= 1 




Thickness of magnetic shell [R„] 


I 


0.3 


-30 


Aloj 


;/ = 


0.05 


Relativistic electron density [cm -3 ] 


fit 


1- 


10 7 


Aloj 


;«r = 


= 0.25 


Relativistic electron energy power-law index 


6 


2 


-3 


A6 = 


0.5 




Temperature of the trapped thermal plasma at stellar surface [K] 


T 


10 4 


-10 7 


Aloj 




= 1 


Electron density of the thermal plasma at stellar surface [cm -3 ] 


"po 


10 s - 


-10" 


Alo£ 




= 0.225 


Free parameters - step 2 


Cross-sectional size of the Cold Torus [R,] 


d a 


1 


-6 


Ad« 


= i 




Temperature of the plasma of the Cold Torus [K] 


T* 


10 3 


-10 5 


Alof 




= 0.25 


Electron density number of the plasma in the Cold Torus [cm -3 ] 


n a 


10 9 - 


-10 12 


Alo£ 




= 0.5 



n r total number density of non-thermal electrons; 
5 spectral index of non-thermal electron energy distribution; 
T Po temperature, at the stellar surface, of the thermal plasma 

trapped in the inner-magnetosphere; 
n Po electron number density, at the stellar surface, of thermal 
plasma; 

d ct cross-sectional diameter of the cold torus surrounding the 
star; 

T ct temperature of the cold torus; 

n a electron number density of the cold torus. 

The range of values and the assumed simulation steps of 
these free parameters are reported in Table|3] 

In TLU04 the Alfven surface has been located numerically 
by solving the equation of balance between kinetic and mag- 
netic energy density. In this work r& is assumed as a free pa- 
rameter, since the wind parameters of CU Vir are still unknown. 
The mass loss rate M, assumed flowing from the whole stellar 
surface, and the plasma density of the wind at the Alfven sur- 
face are thus derived as output parameters for a given r& and 
wind velocity v«,. 

For a detailed discussion of the computational method we 
refer to TLU04. 

4.3. Numerical simulations 

Since the computation is quite time consuming, we split the 
analysis in two steps. 

4.3.1 . Step 1 : fitting the 8.4 GHz light curve 

In this first step we can neglect the presence of the cold torus. 
As pointed out by TLU04, it does not have an appreciable effect 
in the low frequency radio emission (5 and 8.4 GHz) that arises 
from external magnetospheric regions, and we exclude a priori 
its effect in the simulations at 8.4 GHz, reducing the number of 
free parameters (Table [3]- step 1). 

For the choice of the sampling length of the three dimen- 
sional grid we adopt the same criteria as in TLU04, with a nar- 
row spacing (0.08 R*) for the distance from the star center of 
less than 2.3 R„, middle spacing (0.3 R„) between 2.3 and 7 R„, 
and rough spacing (1 R.) greater than 7 R*. 



Table 4. Solutions of the 3-D model: the free parameters able 
to reproduce the observations. 



Symbols 




Solutions 


Units 


r A 




12-17 


[R.] 


S 




2 




da 




5 


[R.] 


l/r A 




0.05 - 0.2 




«r 


oc rjl 


0.88 - 3.54 


[xlO 3 cm -3 ] 


«P0 


- 1 fl5.95 7*0,68 
PO 






Bet 


> 10 595 7- c ° t 68 







We compute the radio emission for the Stokes / and V at 
the phases that characterize the observed light curve at 8.4 GHz 
(see Sect. QJ- The match between simulations and observations 
defines the possible configurations of the magnetosphere. 

4.3.2. Step 2: fitting the spectra 

For any set of parameters modeling the observed Stokes / and 
V at 8.4 GHz we compute the radio emission at 5 and 15 GHz 
to search for the ones fitting the whole observed spectra and 
their polarization degree. 

However we cannot rule out the importance of the cold 
torus on the stellar emission at the frequency of 15 GHz. The 
free parameters that characterize the cold torus and the adopted 
ranges are reported in Table|5]- step 2. 

5. Results 

5.1. Prediction of the simulations 

Comparison between the observations and simulations allow 
us to find the combinations of the model free parameters that 
can account for the observational behavior. The convergence of 
these quantities gives important informations about the physics 
of the external regions of this star. 

The solutions derived from our analysis are reported in 
Table the corresponding observed and computed flux light 
curves are plotted in Fig.|3] 
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Fig. 3. Comparison between CU Vir observed and computed flux density and degree of circular polarization, as a function of 
rotational phase. 



The free parameters 5 and d ct are determined with an un- 
certainty smaller than the simulation steps. We find acceptable 
solutions for r A in the range 12 - 17 R», / in 0.05 - 0.2 r A 
and n r in 0.88xl0 3 - 3.54xl0 3 cm 4 . We note that the rela- 
tion n, oc r A /l holds (Table |4). An important conclusion can 
be drawn: the column density of the relativistic electrons at the 
Alfven radius is uni vocally determined by r A . For the different 
values of the Alfven radius deduced from our simulations the 
column density is in the range 3.2 x 10 14 - 4.6 x 10 14 cm . 

Regarding the physical properties of the thermal plasma 
trapped in the inner magnetosphere, we found that the values 
of the electron number density and temperature that account 
for the observed spectral index are related as 



_ 1Q5.95 71O.68 
"po - iU 1 po ' 



(1) 



while for the plasma in the cold torus 

„ ^ 1 r\5.95 t-0.68 

In terms of free-free absorption coefficients (Dulk 1985 1, 



k v = 9.78 10 



-3 n z 



18.2+ fin T- In v (T < 2 10 5 K) 
24.5 + In T - In v (T > 2 10 5 K). 



The above relations imply that the plasma surrounding the star 
is almost opaque at v =15 GHz, and that the cold torus around 
the star is completely opaque at the same frequency. 

The Alfven radius ta is a free parameter of our model. It is 
defined as the point where the magnetic energy density is equal 
to the kinetic energy density. If the flowing wind has terminal 
velocity Vco and density p w = oth« w , with proton mass mn, 
then: 



B 1 



1 



For a given r^, we obtain the corresponding magnetic field 
strength B(r A ), and consequently the thermal plasma density of 
the wind HwC^a) at the Alfven surface. In Table [5] we list the 
resulting values for « w ('"a) as a function of r& assuming three 



different values of the wind velocity (v„ = 400, 600 and 800 
km s~'). 

The mass loss rate for a spherically symmetric wind is eas- 
ily computed from the continuity equation 

M = 47ir A v 00 mH«w 

In the presence of a magnetic dipole, only the wind flowing 
out from the polar caps can escape from the star and the actual 
mass loss rate M act is a fraction of M, defined as the ratio of the 
area of the polar caps to the whole stellar surface (see TLU04). 
The results for M and M act are shown in Table|5] 



We rejected all the solutions that do not satisfy the condi- 
tion of equality between the thermal plasma energy density and 
the wind ram pressure in the equatorial region of the Alfven 
surface. 

As previously discussed, the filling of the inner magneto- 
sphere by thermal plasma is a consequence of the radiatively 
driven stellar wind, which accumulates material until a steady 
state, along lines of force of the magnetic field, is established, 
i.e. the thermal plasma pressure e p balances the wind ram pres- 
sure pram. This condition at the Alfven point can be written as: 

— Pram 



k B n v (r A )T p (r A ) = m H n w (rA)v , 



(2) 



The stellar rotation implies « p = n Pa r~ l and T p = T Po r, hence e p 
is constant inside the inner magnetosphere, whereas the p ra m is 
defined by the wind velocity Voo and wind density at r A (Babel 
& Montermle 1997 1. Combining Eq.|5]and Eq.[2 we can derive 
the temperature and density of the thermal plasma trapped in 
the inner magnetosphere. The so-derived values of T pa and n po 
are reported in Table|5](last two columns). 

In Fig. 0] the shaded area represents the pairs of log T and 
logn for which the thermal plasma is optically thick at v =15 
GHz (cold torus), while the Eq.[^and Eq.|2] for different values 
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Table 5. Results for wind, mass loss and inner magnetosphere: 
as a function of the Alfven radius r&, the spherical wind mass 
loss rate M, the effective mass loss M act and the number density 
of the electrons of the wind at the Alfven point Hw^a) are re- 
ported. In the last two columns, temperature (r po ) and number 
density (n Po ) at the base of the magnetosphere are also reported. 
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of ta, are represented by the continuous and dashed line respec- 
tively. The intersection, marked with a thick segment, defines 
the ranges of values of T Po and n Po that satisfy the equality con- 
dition and correctly reproduce the observed radio light curves. 

5.2. Dependence on the stellar parameters 

Since the distance and stellar radius of Cu Vir are known with 
good accuracy, the major uncertainty on the assumed param- 
eters is related to the magnetic field strength and geometry. 
In particular, the value of B p reported in the literature ranges 
from 2200 to 3800 Gauss (Landstreet et al. 1977; TLLOO). We 
thus examined whether our choice of B p =3000 Gauss signifi- 
cantly influences the model results. We verified that a variation 
of ±25% in B p results in a variation of +25% in « r , while the 
other parameters remain unaffected. 

Moreover, we tested the sensitivity of the model to values 
of inclination of the rotation axis (z) and obliquity of the mag- 
netic axis (B) different from those adopted in our calculations. 
We minimize the normalized x 2 function for the light curves 
(total flux density and the fraction of circular polarization) at 
8.4 GHz. To restrict the computational time, we limit this anal- 




1'3 



io 3 io 4 io 5 

T [K] 

Fig. 4. log T - log n plane for the magnetospheric plasma. In the 
shaded area the cold torus has Tisghz > 1 ■ For the inner mag- 
netosphere, the thermal pressure along the dashed lines equals 
the wind ram pressure, as computed for the possible values of 
^a. The thick segment of the continuous line (Eq.^ represents 
the loci where the T Po and n Po of our model reproduce the ob- 
servations. 



ysis to the set of model parameters associated with the lowest 
value of the Alfven radius (rA = 12 R„). The maps of^- 2 are 
displayed in Fig.|5](top panels). Dark areas correspond to pairs 
of i and B with lower x 2 , for which we find the best match be- 
tween simulations and observations. We note that the x 1 from 
the total flux fit is not well correlated to the x 1 relative to the 
fit of the circular polarization. The product of these two ma- 
trices, Fig. |5] (bottom panel), gives the combination of i and B 
that better reproduces both total flux and circular polarization. 
We note that the most probable values of i and B are restricted 
in a well defined range of values, respectively: i = 41° - 52° 
and B = 68° - 77°. We note that the values of i and B from the 
literature are well inside the previous intervals. 

5.3. Limits of the model 

While the modulation of the total flux density light curve ob- 
served at 5 GHz is well reproduced by our simulations, we 
are not able to match the two peaks observed close to the 
phases 0.40 and 0.85 of the circularly polarized flux light curve. 
These are also the phases of the coherent emission observed at 
1.4 GHz (Fig.n bottom panel), a coincidence suggesting that 
the two peaks are residuals of the Electron Cyclotron Maser 
Emission (ECME) mechanism (TLLOO). As the ECME mech- 
anism is not considered in our model, which computes only in- 
coherent gyrosynchrotron, it cannot reproduce these observed 
peaks. 

At 15 GHz, a stronger and broader polarization is expected 
in the range O = 0.9 - 1.3 (Fig. [3}. This discrepancy could 
be ascribed to the simple dipolar configuration adopted for the 
magnetic field. In the magnetospheric regions close to the star, 
where the bulk of 15 GHz radiation is emitted, the high or- 
der components of the magnetic field of CU Vir (Hatzes [l9971 
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Fig. 6. Simulated maps for stellar orientations associated with the radio light curve extrema. On the left side the corresponding 
radio maps are displayed; the theoretical bightness spatial distribution is averaged over half the rotational period and convolved 
with the 1 mas beam of the VLBI array. The contours are 10%, 20%, 30%, 50%, 70%, 90% of the peak brightness. 




30 40 50 60 70 30 40 50 60 70 

i i 




Fig. 5. Top panels: Maps, as a functions of i and /?, of the nor- 
malized^ 2 obtained for the light curves at 8.4 GHz, total flux 
(left) and fraction of circular polarization (right). The^- 2 min- 
ima are respectively w 2 and » 1.6. Bottom panel: product of 
the two top panels. 



Kushinig et al. [T999l TLL00; Glagolevskij & Gerth 2002 1 can- 
not be neglected. A simple dipole provides a symmetric total 
flux light curve and an enhancement of the fraction of polar- 
ization when the line of sight is close to the magnetic pole, 
where the magnetic lines are almost parallel. In the presence 
of multipolar components the radio flux light curves would not 
be symmetric and the fraction of the polarization would be re- 
duced. 

6. Discussion 

In the following we will analyze physical information about 
the stellar magnetosphere that can be inferred from the results 
of our simulations. 

6.1. The acceleration process 

The convergence of the value of the spectral index of the emit- 
ting electron distribution (N(E) oc E~ s ) toward 5-2 indicates 
that the acceleration process occurs in the current sheets just 
outside the Alfven radius. This value agrees with the one de- 
rived in the case of solar flares by Zharkova & Gordovskyy 
( 2005 ), who found a spectral index 5-2 for the non-thermal 
electrons accelerated in the neutral reconnecting current sheets 
of solar magnetic loops. In addition, from the ratio between n r 
(Table |4) and the density of the electrons of the wind at the 
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Alfven point n w (r A ) (Table [3Jl, we can estimate that the effi- 
ciency of the acceleration process for non-thermal electrons is 
2 x 10~ 3±1 . The efficency inferred by our analysis is consis- 
tent with the acceleration efficiency in large stationary unstable 
current sheets associated with coronal loops during solar flares 
(Vlahos et al. l2004l 

6.2. The radio emitting region 

We investigated the possibility of discrimining among the pre- 
vious possible values of the Alfven radius by using high an- 
gular resolution radio observations in the future. In Fig.[6]the 
brightness spatial distributions at 8.4 GHz obtained for the ex- 
treme values of r A (respectively 12 and 17 R t ) are shown. The 
simulated radio maps have been calculated for three rotational 
phases: the stellar orientations that show the north and south 
poles, associated with the two maxima and the orientation that 
gives a null effective magnetic field, associated with the mini- 
mum of the radio light curve. In the right panel of Fig. [6] the 
different extension of the radio source for the two cases an- 
alyzed is clearly visible. However such a difference is not so 
obvious when looking at the radio emission arising from re- 
gions near the stellar surface, where the highest fraction of the 
magnetospheric radio emission originates. 

For the Very Long Baseline Interferometry (VLBI) tech- 
nique, that can be used to map CU Vir at milli arcsecond scales, 
we derived two "synthetic" maps at VLBI resolution, shown on 
the left side of Fig.|6j obtained as the average of hypothetical 
8 hours observation around the stellar phase coinciding with 
the null effective magnetic field, convolved with a Gaussian 1 
mas beam. It is clear that the present capability of VLBI obser- 
vations does not allow us to discriminate between the extreme 
values of the Alfven radius of CU Vir predicted by our model. 



6.3. The thermal plasma 

The physical condition in the inner magnetosphere play an im- 
portant role in the emerging radio flux because of the free-free 
absorption by the plasma, which affects the modulation of the 
radio light curves and the spectrum of the source. The compar- 
ison between our model prediction and multifrequency obser- 
vations reveals the presence of a torus of material colder and 
denser than the trapped thermal plasma. This is necessary to 
reproduce the correct scale of the total flux density in the 15 
GHz light curve of CU Vir. The electron number density (« ct ) 
and temperature (T ct ) of the plasma of the cold torus, are within 
the shaded area of the plane log T-log n shown in Fig.|4] where 
the optical depth of the cold torus is higher than 1 at 15 GHz. 
Higher frequency measurements (for example 22 and 43 GHz) 
could characterize the spatial distribution of the thermal plasma 
in the deep magnetospheric regions, near the stellar surface. 

The temperature and density of the thermal plasma trapped 
in the inner magnetosphere are assigned for a given Alfven ra- 
dius and wind velocity (Table[5)- As a by-product of our anal- 
ysis, we estimate the X-ray emission from the thermal plasma 
trapped in the stellar inner magnetosphere of CU Vir as a func- 
tion of r A and v M . Following the same procedure described in 
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Fig. 7. Predicted X-ray fluxes of the inner magnetosphere of 
CU Vir for the solutions of our 3-D model. 



TLU04, for each element of the 3D cubic grid that samples 
the magnetosphere, we computed the thermal bremsstrahlung 
emission coefficients and the emerging power integrated be- 
tween 0.1 and 10 keV, in the typical energy range of X-ray 
telescopes like Chandra and XMM. The resulting X-ray fluxes 
at the Earth are displayed in Fig. [71 

The stellar X-ray emission is a function of the size of the 
Alfven surface (Fig. [7J| and, in particular, it decreases as the 
Alfven radius increases. This is a consequence of the lower 
temperature and density of the plasma in the inner magneto- 
sphere for increasing r A . For the same reason, lower wind ve- 
locities cause a lower X-ray emission. 

On the basis of the values derived from our calculations, 
we conclude that the X-ray emission from CU Vir is below 
the detection limit of the ROSAT all sky survey (the average 
background at the position of CU Vir in the RASS is higher 
than 10~ 12 erg ctrT 2 s _1 ). Nevertheless, for exposure times of a 
few tens of kilo-seconds, Chandra and XMM can reach an X- 
ray detection limit close to 10~ 15 erg cirr 2 s , which is lower 
than the X-ray fluxes predicted from several model solutions. 

Within our model, measurement of the X-ray flux of CU Vir 
could place an important constraint on the physical condition 
in the magnetosphere. 



7. Conclusion 

The radio magnetosphere of CU Vir has been extensively stud- 
ied using of the first multifrequency VLA observations, cov- 
ering the entire rotational period, for both the total flux den- 
sity (Stokes /) and the circular polarization (Stokes V). We have 
successfully modeled the observed flux modulation and the cir- 
cular polarization light curves using the 3D model developed 
by us (TLU04), modified to resolve the circular polarization 
intensity. 

We pointed out the importance of multi-frequency and dual 
polarization observations throughout the rotational period to 
gain insight into the physical conditions of the circumstellar 
environments surrounding MCP stars. 
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Our results are: 

- We confirm the physical scenario proposed by Andre et al. 
fT988l and Linsky et al. {1992) and the MCWS model pro- 
posed by Babel & Montermle J1997I : 

- The size of the CU Vir inner magnetosphere has been con- 
strained in the range: r A = 12 - 17 R«; 

- we evaluated the mass loss rate of CU Vir, which up to now 
was not determined for this star, to be about 10~ 12 Mq yr _I ; 

- the temperature of the thermal electrons trapped in the in- 
ner magnetosphere can reach 10 5 - 10 6 K and can give de- 
tectable X-ray emission; 

- the presence of an accumulation of cold and dense material 
around the star has been shown and limits for its size and 
physical conditions (temperature and density) are given; 

- the acceleration in the current sheets produces a non- 
thermal electron population with a hard energetic spectrum 
(N(E) oc E~ 5 with 6 = 2), and the acceleration process has 
an efficiency of about 1CT 3 . 

The limits of our model to reproduce all the observed fea- 
tures, for example the two peaks at phases « 0.4 and « 0.85 
observed in the 5 GHz light curve for the circular polarization 
and the shape of the light curves at 15 GHz, have been inter- 
preted respectively as due to the residual of the ECME mech- 
anism (responsible for the features previously observed at 1 .4 
GHz, TLL00) and as evidence of a multipolar component in 
the magnetic field of CU Vir. 

Acknowledgements. We thank the referee for his/her constructive crit- 
icism which enabled us to improve this paper. 
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Appendix A: Calculation of circular polarization 

The circularly polarized intensity may be obtained using the 



(A.l) 



relation, given by Ramaty ( 1969 ): 



V = 2 



1 + Qg + 2 1 + dgj 1 



where I ± are solutions of the radiative transfer equation for the 
ordinary (sign "+") and extraordinary (sign "-") modes (Klein 
& Trotter 1984): 



d I ± \cosif/ ± 

ds llg ± 2 



| COS tff ± 

ne± 2 



(j] ± - k ± I ± ) 



(A.2) 



with tj ± and k ± gyrosynchrotron emission and absorption co- 
efficients, ng ± and ag ± refractive index and polarization coeffi- 
cient for the two modes respectively. if/ ± is the angle between 
the wave vector and the group velocity, given by the equation: 



iff ± = arctan 



dn 



ng ± d6j 



with Q the angle between the wave vector and the magnetic 
field. The equation |A. 21 was numerically solved, along a path 
parallel to the line of sight, using the method developed by us 
inTLU04. 

In accordance with classical physics, for propagation paral- 
lel to the magnetic field lines (6 = and a g± — ±1), the ordinary 
mode is characterized by the electric field vector rotating in the 
sense of the electrons, counterclockwise in the case of a wave 
approaching the observer, therefore the polarization sense for 
the ordinary mode is LCP; conversely the extraordinary mode 
has a polarization sense RCP. 

VLA measurements of the wave circular polarization state 
are in accordance with the IAU and IEEE orientation/sign con- 
vention, unlike the classical physics usage. Therefore to com- 
pare our simulations with the VLA observations we have to 
change the sign of the result of eauation lA.il 



